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Abstract 

Quantum amplitudes for s = 1 Maxwell fields and for s — 2 lin- 
earised gravitational-wave perturbations of a spherically-symmetric Ein- 
stein/massless scalar background, describing gravitational collapse to a 
black hole, are treated by analogy with the previous treatment of s = 
scalar-field perturbations of gravitational collapse at late times. Both the 
spin-1 and the spin-2 perturbations split into parts with odd and even 
parity. Their detailed angular behaviour is analysed, as well as their 
behaviour under infinitesimal coordinate transformations and their lin- 
earised field equations. In general, we work in the Regge- Wheeler gauge, 
except that, at a certain point, it becomes necessary to make a gauge 
transformation to an asymptotically-flat gauge, such that the metric per- 
turbations have the expected fall-off behaviour at large radii. In both the 
s = 1 and s — 2 cases, we isolate suitable 'coordinate' variables which can 
be taken as boundary data on a flnal space-like hypersurface Ef • (For 
simplicity of exposition, we take the data on the initial surface E/ to be 
exactly spherically-symmetric.) The (large) Lorentzian proper-time inter- 
val between E_r and Ef , measured at spatial inflnity, is denoted by T . 
We then consider the classical boundary- value problem and calculate the 
second- variation classical Lorentzian action S^^^^^ , on the assumption that 
the time interval T has been rotated into the complex: T — > \T\ exp(— i9), 
for < 9 < tt/2 . This complexified classical boundary- value problem is 
expected to be well-posed, in contrast to the boundary-value problem in 
the Lorentzian-signature case (9 = 0), which is badly posed, since it refers 
to hyperbolic or wave-like field equations. Following Feynman, we recover 
the Lorentzian quantum amplitude by taking the limit as 9 — > 0+ of the 
semi-classical amplitude exp^iS'^j^^g^^ . The boundary data for s = 1 in- 
volve the (Maxwell) magnetic field, while the data for s — 2 involve the 
magnetic part of the Weyl curvature tensor. These relations are also in- 
vestigated, using 2-component spinor language, in terms of the Maxwell 
field strength (pAB = 0{ab) and the Weyl spinor 'I'ascd = '^(abcd) ■ 
The magnetic boundary conditions are related to each other and to the 
natural s — ^ boundary conditions by supersymmetry. 
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1 Introduction 



This paper describes part of a project concerned with the calculation of quan- 
tum amplitudes (not just probabilities) associated with quantum fields, includ- 
ing gravity itself, in the case that strong gravitational fields may be present. 
The most obvious example - the original motivation for this work concerns 
quantum radiation associated with gravitational collapse to a black hole [1-f 1]. 
But the framework adopted here is more general, and certainly does not depend 
on whether there is a classical Lorentzian-signature collapse to a black hole. It 
includes the case of local collapse which is not sufficient to lead to (Lorentzian) 
curvature singularities, and also quantum processes in cosmology, where, for ex- 
ample, anisotropics in the Cosmic Microwave Background Radiation (CMBR) 
can be computed, and depend crucially on the underlying Lagrangian for gravity 
and matter [12]. 

To exemplify the underlying ideas, we consider the case of local collapse 
(whether or not to a black hole). Thus, the gravitational field is taken to be 
asymptotically flat. For simplicity, consider Einstein gravity coupled minimally 
to a massless scalar field . In classical gravitation, we are used to describing 
this by means of a Cauchy problem, giving evolution to the future (say) of an 
initial space-like hypersurface S , which extends to spatial infinity. We write 
g^i, {HtV = 0,1,2,3) for the components of the 4-dimensional metric, and 
then denote by hij = gij = 1, 2, 3) the components of the intrinsic spatial 
metric on S in the case that S is given by the condition x° = const . Cauchy 
data would, loosely speaking, consist of hij , (j) and their corresponding normal 
derivatives on <S. By contrast, in quantum theory, one typically asks for the 
amplitude to go from an initial configuration such as {hij on an initial 
hypersurface E/ , to a final configuration {hij, (f>)p on a final hypersurface Tip . 
The problem of finding the quantum amplitude should (naively) be completely 
posed, once one has also specified the (Lorentzian) proper-time interval between 
the surfaces S/ and Sj- , as measured near spatial infinity. 

Much of the 'non-intuitive' nature of quantum mechanics can be traced to 
the 'boundary- value' nature of such a quantum amplitude [13], as compared 
with the familiar classical initial-value problem. A crucial mathematical aspect 
of this difference, responsible for a good part of the 'non-intuition', is that the 
classical version of the problem of calculating a quantum amplitude, as posed 
above, would involve solving the classical field equations (typically hyperbolic), 
subject to the given boundary data {hij ,4')i,f on the hypersurfaces ,Y,p , 
separated near spatial infinity by a Lorentzian time interval T . As is well known, 
a boundary-value problem for a hyperbolic equation is typically not well posed. 
For typical boundary data, a classical solution will not exist 14,15]; or, if it 
does exist, it will be non- unique. The straightforward cure for this ill, due to 
Feynman [13], is of course to rotate the Lorentzian time-interval T into the 
complex: T \T\ exp(- iO) , with < 6* < 7r/2 . 

A simplified classical boundary-value example, showing this behaviour, is 
described in [16]. When this example is posed originally (and badly) in 2- 
dimensional Minkowski space-time, one considers a scalar field 0(t , x) , obeying 
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the wave equation 



--g^ + g^ = ; 0<t<T , -oo<a;<+(X) , (1.1) 

on the assumption that (j) decays rapidly as |a;|— >oo. A simple choice of 
Dirichlet boundary data is to take 

(j){t=0,x) = , (j){t = T ,x) = (j)x(x) . (1.2) 

The time-interval T at spatial infinity is then, as above, rotated into the 
lower complex half-plane: 

T \T\ exp(-i6i) ; < 6» < 7r/2 . (1.3) 

For convenience, we define, for a given fixed 6 (0 < ^ < 7r/2) , the 'rotated- 
time' coordinate 

y = texp{-i9) . (1.4) 
In terms of the new coordinates {y ,x) , the wave equation (1.1) reads 

and the boundary conditions (1.2) become 

cj>{y = Q,x) = , (l){y = T^\x) = M^) ■ (1-6) 

Here, the extreme case 6 = tt/2 corresponds to the Riemannian (Euclidean) 
sector, and to a well-posed real elliptic Dirichlet boundary-value problem for 
the Laplace equation. The potential (p{y , x) is thus required to be a complex 
solution of Eqs.(1.5,6), which, for < 9 < tt/2 , describe a strongly elliptic 
partial differential equation in the sense of [17]. The property of strong elliptic- 
ity guarantees existence and uniqueness in this linear example. As verified in 
[16], however, the 'classical solution' becomes singular in the Lorentzian limit 
0^0+ . 

In our coupled non-linear gravitational/scalar-field example, the extreme 

case 9 = tt/2 would correspond to a purely Euclidean time-interval |r| , and 
classically one would then be solving the field equations for Riemannian gravity 
with a scalar field (p . Since these field equations are 'elliptic modulo gauge' - see 
[17] - one would expect to have a well-posed classical boundary- value problem, 
with existence and uniqueness. The intermediate case < 9 < tt/2 requires 
the interval T and any classical solution to involve the complex numbers non- 
trivially. If the problem turns out to be strongly elliptic, up to gauge, then the 
complex case < 9 < tt/2 would again be expected to have the good existence 
and uniqueness properties of the real elliptic case. 

In practice, in the black-hole evaporation problem or (say) in cosmological 
examples, one typically treats the case in which both the gravitational and the 
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scalar initial data arc close to spherical symmetry. Hence, as a leading ap- 
proximation, one begins by studying the spherically-symmetric Einstein/scalar 
system. This was treated in [18] for Lorentzian signature and is outlined in [19] 
for Ricmannian signature. In the Riemannian case, the metric is taken (without 
loss of generality) in the form 

ds^ = e'' dT^ + dr'' + {de^ + sin'' 9 dcj)'') , (1.7) 

where 

b = 6(r,r) , a = a(T,r) , (1.8) 

and the scalar field is taken of the form (f){T,r) [19]. The scalar field equation 
reads: 

+ e''-'^0" + i (a -&)</) + r-ie^-'^(l + e'')<^' = , (1.9) 

where (' ) denotes d{ )/dT and ( )' denotes d{ )/dr. Together with Eq.(1.9), a 
slightly redundant set of gravitational field equations is given by: 

a' = - Awr {e''"'' - <j>'^ ) + r-'^ {1 - e"") , (1.10) 
b' = -47rr(e''-''</)2 -</)'^) -r-i(l-e") , (1.11) 
a = Snr^fl)' , (1.12) 

a+e^-°6" + i(d-6)d-r-^e''-"(l-e")(6' + 2r-i) = Stt (^Ve^"" (A'^). (1.13) 

The metric and the classical field equations in Lorentzian signature [18] can 
be derived from the above by the formal replacement 

t = re-'^ , (1.14) 

where 9 = it/ 2 is independent of 4-dimensional position. Similarly for complex 

metrics with suitable behaviour at infinity, with < 6 < 7r/2 . 

Even in the spherically-symmetric case, very little is known rigorously about 
existence and uniqueness for the Riemannian (or complex) boundary- value prob- 
lem. For this case, numerical investigation of the weak-field Riemannian boundary- 
value problem was begun in [19], and has recently been extended towards the 
strong-field region [20]. For weak scalar boundary data, global quantities such 
as the mass M and Euclidean action / appear to scale quadratically, in accor- 
dance with analytic weak-field estimates [20] . In the limit of strong-field scalar 
boundary data, it may be that a typical pattern will emerge numerically for the 
general 'shape' of the classical Riemannian gravitational and scalar fields. In 
that case, it might be possible to find analytic approximations for the strong- 
field limit (quite different from those valid in the weak- field case), which could 
provide further analytical insight into the solutions of the coupled Riemannian 
Einstein/scalar boundary-value problem. In particular, it would be extremely 
valuable to have strong-field approximations which were valid into the complex 
region, with < 9 < 7r/2. One might conjecture that, as one approaches 
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the Lorentzian limit 0^0+, for very strong spherically-symmetric bound- 
ary data, the solutions (albeit complex) correspond to classical Einstein/scalar 
solutions which form a singularity, surrounded by a black hole. 

In the case of quantum amplitudes for Lagrangians with Einstein grav- 
ity coupled to matter, one can consider anisotropic boundary data posed in 
the 'field language' of this paper, by taking (for the present Einstein/scalar 
case) non-sphcrically-symmetric boundary data {hij , 4>)i,f ■ Then, at least in 
the asymptotically-flat case with time-interval T at spatial infinity, one is in- 
evitably led to consider the complexified boundary-value problem, with T — > 
|T| exp(— 20) . but with imchangcd data [hij ,<p)i^F on the other boimdaries. 
This corresponds (by a slight re-definition oi 6) to the procedure adopted in 
our model 2-dimensional boundary-value problem of Eqs. (1.1-6). Even for fairly 
small 9 , solution of this boundary- value problem is expected to smooth out vari- 
ations or oscillations of the boundary data, when one moves into the interior by 
a few multiples of the relevant wavelength. If the problem is genuinely strongly 
elliptic, up to gauge, then one will be able to extend the classical solution ana- 
lytically into the complex. 

Strictly, in order that quantum amplitudes should be meaningful for any 
Einstein- gravity /matter Lagrangian imder consideration, one should work only 
with theories invariant under local supersymmetry - that is, with supergravity 
models or supergravity coupled to supermatter [21,22]. Thus, for example, the 
bosonic Einstein/massless-scalar model above should be replaced by the simplest 
locally-supersymmetric theory which contains it [21]. This N = 1 supergrav- 
ity/supermatter model contains a complex scalar field (f) , with a massless spin-^ 
partner; the graviton acquires a spin-| gravitino partner. Generally, for 'Ric- 
mannian' boundary data, one expects that the resulting 'Euclidean' quantum 
amplitude has the semi-classical form 

Amp ~ {An + HA^ + h^A2 + ...) exp(- /dass/^) , (1-15) 

asymptotically in the limit that (/class /^) — > . Here, /class is the classical 'Eu- 
clidean action' of a Riemannian solution of the coupled Einstein and bosonic- 

matter classical field equations, subject to suitable boundary conditions. In the 
complex regime of this paper, we shall use the expressions / and — iS inter- 
changeably, where S denotes the 'Lorentzian action'. For simplicity, we assume 
that there is a unique classical solution, up to gauge, coordinate and local super- 
symmetry transformations. But it is quite feasible, in certain theories and for 
certain boundary data, to have instead (say) a complex-conjugate pair of clas- 
sical solutions [23]. The classical action /class and loop terms Aq , Ai , A2 , . . . 
depend in principle on the boundary data. In the case of supermatter coupled 
to A'' = 1 supergravity, each of /dass ,Ao,Ai,A2,... will also obey differen- 
tial constraints connected with the local coordinate and local supersymmetry 
invariance of the theory, and with any other local invariances such as gauge 
invariance (if appropriate) [15,24]. 

In particular, in the locally-supersymmetric case, the semi-classical expan- 
sion (1.15) may become extremely simple [15,25,26]. For example, for N = 1 
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supergravity, for purely bosonic (Einstein) boundary data, one has [15]: 

Amp ~ Ao exp(-/ciass/?i) , (1.16) 

where, further, the one-loop factor Aq is in fact a constant. When the boundary 
data are allowed to include both bosonic and fermionic parts (suitably posed), 
one expects that a classical solution of the coupled bosonic/fermionic field equa- 
tions will still exist. In this case, the expression /class in Eq.(1.16) denotes the 
full classical action, including both bosonic and fermionic contributions. The 
fermionic contributions (naturally) also depend on the boundary data, and, as is 
standard in the holomorphic representation used here for fermions [27-29], live 
in a Grassmann algebra over the complex numbers. Related properties hold for 
N = 1 supergravity coupled to gauge- invariant supermatter [22,25]. Whether 
or not the supermatter is also invariant under a gauge group, there will be anal- 
ogous consequences for the semi-classical expansion (1.15) of the quantum am- 
plitude. In particular, one expects finite loop terms Aq ,Ai ,A2 . . [15,25,26], 
but typically not the maximal simplicity of the pure supergravity amplitude 
(1.16). 

In the case (1.16) of pure supergravity, the classical action is all that is needed 
for the quantum computation. A corresponding situation arises with ultra-high- 
energy collisions, whether between black holes [30], in particle scattering [31], 
or in string theory [32]. 

To fix one's physical intuition, one can assume that, near the initial surface 
E/ , the gravitational and scalar fields are approximately spherically symmetric 
and vary extremely slowly with time, corresponding to diffuse bosonic matter 
near E/ . The final hypersurface E f should preferably be taken at a sufiiciently 
late time T that all the quantum radiation due to the evaporation of the black 
hole will by then have been emitted. The final data {hij , ^) f for gravity and the 
scalar field, together, if a Maxwell field is included, with the spatial components 
{Ai)F of the vector potential, are taken to have small anisotropic parts - this 
corresponds, in 'particle language', to a choice of final particle state. The re- 
sulting 'weak-field' quantum amplitude, to be found below, can be described in 
terms of products of zero- or one-particle harmonic-oscillator eigenstates. One 
could, of course, also consider final data which deviate strongly from spher- 
ical symmetry. Their quantum amplitudes will still be roughly proportional 
to the (complex) quantity lim0^o+ cxp(— /dass) • For the weak perturbations 
(above) away from spherical symmetry, /dass is nearly quadratic, but for strong 
perturbations, /dass will be very non-linear, and the probabilities of such final 
configurations will be microscopic. 

In considering the classical boundary- value problem, even though spin-1 and 
fermionic fields may also be involved, we shall for simplicity first consider the 
fields and (p . As above, we consider the Riemannian (or complex-rotated 
Riemannian) classical boundary- value problem, given boundary data {hij ,4>)i 
and (hij ,(f>)F , where the initial and final boundary hypersurfaces E/ and Ei? 
are separated at spatial infinity by a complex (Riemannian-time) interval of 
the form |T| exp(— i9) . Since all fields are to be regarded as perturbations of 
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a 'background' spherically-symmetric configuration, we are assuming that the 
classical solutions {g^i, , (j)) of the coupled Einstein/scalar field equations may be 
decomposed into a 'background' spherically-symmetric part (7^,^ ,$), together 
with a 'small' perturbative part. The linearised pcrturbativc fields, whether 
spin-0 scalar [11,33,34], spin-1 or spin-2 (this paper), spin-^ [35] or spin-| (in 
progress [36]), can be expanded out in the appropriate spin- weighted spherical 
harmonics [37-39]. 

Typically, the perturbative scalar-field configuration (say) (pp^ , given on the 
late-time surface Sj? , will involve an enormous number of modes, both angular 
and radial, but with a minute coefficient for each mode. That is. the given t/)^^ 
may contain extremely detailed angular structure, and also be spread over a 
considerable radius from the centre of spherical symmetry of the background 
(7^,^,$), again with detailed radial structure. Similar comments should ap- 
ply to the perturbative part hijp of the spatial gravitational field on T,p , and 
to the final spatial spin-1 potential Aip , if appropriate. As a result, in the 
(complexified) nearly-Lorentzian regime, one will have (classically) radiation of 
various spins, typically with wavelengths much shorter than the characteristic 
length- or time-scale corresponding to the Schwarzschild mass M ; this radia- 
tion will propagate approximately by geometrical optics. In turn, the effective 
energy-momentum tensor T^^ due to this radiation will, on the average, be 
nearly spherically symmetric, and will indeed have the form appropriate to a 
radially-outgoing null fluid [40,41]. The classical 'space-time' metric resulting 
from such a null-fluid effective T^;^ is precisely of the Vaidya type [41]. This 
resembles the Schwarzschild geometry, except that the role of the Schwarzschild 
mass M is taken by a mass function m{t , r) , which varies extremely slowly 
with respect both to t and to r in the space-time region containing outgoing 
radiation. The perturbative fields, such as i/)'^' , propagate adiabatically in this 
classical solution. 

In Sec. 2 we shall describe, in tensor language, the boundary data on the 
final surface Ef which are natural for the s = 1 Maxwell and for the s = 2 
graviton cases. These are, respectively, the Maxwell magnetic field Bi and the 
magnetic part Hik of the Weyl tensor. In Sec. 3 we rephrase the s = 1 and 
s = 2 problems in terms of 2-component spinors; the Maxwell field strength 
is determined by the symmetric spinor = (ji^AB) and the Weyl curvature 
by = "^(ABCD) ■ In both cases, the natural boundary data involve 

a 'projection' on symmetric spinors; for s = 1 and 2 , the data are related 
by supersymmetry. In Sec. 4, for the Maxwell field, we begin the process of 
decomposing the classical problem in terms of odd- and even-parity harmonics, 
following Regge and Wheeler. A similar procedure is carried in Sec. 5 for the 
odd-parity gravitational perturbations. Returning to the Maxwell case in Sec. 6, 
the (magnetic) boundary conditions on are described, together with the 
classical action S^^^^ as a functional of the final boundary data. Sees. 7,8 and 9 
are concerned with gravitational perturbations. In Sec. 7, the odd-parity s = 2 
problem is treated roughly by analogy with the Maxwell problem of Sec. 6. Even- 
parity gravitational perturbations are treated in Sees. 8 and 9. A preliminary 
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analysis in Sec. 8 leads to the classical action functional S^il^^ in Sec. 9. The 
Conclusion is in Sec. 10. 

2 Boundary data for the Maxwell field and for 
gravity 

The Maxwell contribution to the total Lorentzian action S is 

gEM = - -L f d^x{-g)-^F^,F'^-^ , (2.1) 

iDTT 

where F^j/ = Pl/jtv] is the Maxwell field strength, while the space-time metric 
Qfj^i, is assumed here to have Lorentzian signature, with g = det(5^y) < . 
The resulting classical Maxwell field equations are 

V/.F^"' = . (2.2) 

The further condition that F^^^, be derivable from a vector potential , as 

Ff,^ = V^A^ - V^Af, = d^A^ - d^A^ , (2.3) 

may equivalently be written in the form of the dual field equations 

y^{*F^^-) = , (2.4) 

where 

*F^.^ = l^^rFpa (2.5) 

is the dual field strength [42,43]. (The Poincare lemma [44] may be applied, 
since we are working within a manifold M. which may be regarded as a slice of 
, with boundary dM. consisting of two hypersurfaces.) The action (2.1) 
is invariant under Maxwell gauge transformations 

A^ ^ A^+ d^A , (2.6) 

in the interior, where A{x) is a function of position. 

As in [33,34,45] for scalar (spin-0) perturbations of spherically symmet- 
ric Einstein/massless-scalar gravitational collapse, and as in the treatment of 
spin-2 (graviton) perturbations below, we shall need the classical action Sdass , 
namely the action S evaluated at a classical solution of the appropriate (slightly 
complexified) boundary-value problem, as a functional of the boundary data. 
From this, one obtains the semi-classical quantum amplitude, proportional to 
exp(iS'ciass) ) and hence by a limiting procedure the Lorentzian quantum ampli- 
tude. In the present (spin-1) Maxwell case, the classical action 5^^^ resides 
solely on the boundary dM , which consists of the initial space-like hypersurface 
Sj and final hypersurface Sjr . There will be no contribution from any large 
cylinder of radius i?oo — ^ oo , provided that we impose the physically reasonable 
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restriction that, as r ^ cxd , the potential should die off faster than , and 
that the field strength F^i, should die off faster than . That is, we impose 
reasonable fall-off conditions at large r on field configurations, such that the 
action S should be finite. (Compare the usual fall-off conditions for instantons 
in Euclidean Yang-Mills theory [44,48,49].) For the above class of Maxwell field 
configurations, the boundary form of the classical Maxwell action is 

^c'^a^s = - d'xhin.A^F^'^ . (2.7) 

Here, as in Sec.l, hij = gij = 1,2,3) gives the intrinsic Riemannian 3- 
metric on the boundary hypersurface S/ or Ej- , and we write h = det(/iy ) > 
0. Further, n'* denotes the (Lorentzian) unit future-directed timelike vector, 
normal to the space- like hypersurface E/ or E^? . 

Given the 3-1-1 split of the 4-metric 5^1/ at each boundary, due to the 
ability to project vectors and tensors normally using n'^ and tangentially using 
the projector [48] 

h^iu = 9nv + n^n^ , (2.8) 

at the boundary, one can project the potential and field strength F^i, into 
'normal' and 'spatial' parts on E/ and Ei;' . In particular, one defines the den- 
sitised electric field vector on the boundary: 

r = - h^E' , (2.9) 

where 

E" = n^F'"' (2.10) 

obeys rii^E'^ = . Further, in a Hamiltonian formulation [49], when one regards 
the spatial components Ai of the vector potential as 'coordinates', then the 
canonical momentum tt' , automatically a vector density, is given by 

Note that the normal component At = — </? ) where ip is the Maxwell scalar 
potential, is gauge-dependent, but that </? does not need to be specified on the 
spacclike boundaries E/ and E7? , and is indeed allowed to vary freely there' and 
throughout the space-time. Its conjugate momentum therefore vanishes. In the 
gravitational case, analogous properties hold for the lapse function A'' and the 
shift vector A'' [15,49]. 

As described in [49-51], it is natural in specifying a classical boundary- value 
problem for the Maxwell field, with data given on the space-like boundaries E/ 
and Tip and at spatial infinity (with Lorentzian proper-time separation T), to 
fix the spatial magnetic field components, described in densitised form by 

= ^e'^'^Fjk , (2.12) 
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on E/ and E p . The cannot be specified freely on the boundary, but are 
further subject to the (hnear) restriction 

diB^ = . (2.13) 

These components are gauge- invariant, and therefore physically measurable, in 
contrast to those of the spatially-projected vector potential Ai . We shall regard 
the space of such B^{x) , on E/ or E^ , as the 'coordinates' for Maxwell theory. 
From the space-time Maxwell equations (2.2), one also deduces the constraint 

di£' = . (2.14) 

Turning to s = 2 (graviton) perturbations of a spherically-symmetric back- 
ground, we describe the boundary conditions found below to be appropriate 
both for odd- and even-parity vacuum s = 2 perturbations. The most suit- 
able s = 2 boundary data involve prescribing the magnetic part of the Weyl 
curvature tensor CapjS [42,43,52,53] on E/ and on T,p . For simplicity, as 
above, we are taking the gravitational initial data on S/ to be exactly spher- 
ically symmetric ('no incoming gravitons'). Of course, in a large part of the 
space-time, one is nearly in vacuo, the Ricci tensor then obeying i?^/? — , 
whence Ca^-yS — RajS-yS , the Riemann tensor. More generally, the Weyl tensor 
is defined by 

Ra/SiS = Cal3fS + ga['yRs]l3 " 9l3['yRs]a " R 9a[-i9&](3 , (2-15) 

where R = g"'^ Ra0 gives the Ricci scalar, and where square brackets denote 
anti-symmetrisation. 

The algebraic symmetries of the Weyl tensor at a point are summarised by 



These imply that Cap^s has 10 algebraically-independent components at each 
point. At a bounding space-like hypersurface, such as T,p , one can, by analogy 
with the Maxwell case, apply a 3 4- 1 decomposition to the Weyl tensor Cap-ys , 
which splits into two symmetric trace-free spatial tensors, the electric part En- 
and the magnetic part Hn- of the Weyl tensor [52,53]. Thus, the 10 space-time 
components of C^p-ys have been decomposed into the 5 spatial components of 
and 5 more of Hi]^ . (Correspondingly, in Maxwell theory above, the 6 non- 
trivial components of the field strength F^^ became the 3 of Ei plus the 3 of 
Bi.) 

For convenience of exposition, consider an 'adapted' coordinate system (.t" ^x^ , x^) 
in a neighbourhood of Ef , such that Ef lies at x° = , and such that -nP = 1 
at all points of Ef • The spatial 3-metric is, as usual, denoted by hij , and we 
again write h = det(/iy ) . The electric part of the Weyl tensor is defined in 
4-dimensional language to be 

Ea-y = Cc^ffjsnf^n^ ■ (2.17) 
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In an adapted coordinate system, this corresponds to the 'spatial' equation 

Eik = Cioko ■ (2-18) 
The magnetic part of the Weyl tensor is defined to be 

= IvarCp^isnPn^ , (2.19) 

where 

VafftS = ViaP-fS] = {-g)^ ^a0j5 (2.20) 

is the alternating tensor, with g = dct(gp,y) , and where eap-^s = ^[a/j^ij] is the 
alternating symbol, normalised such that 60123 = 1 . In an adapted coordinate 
system, one finds that 

Hik = — - hin e"^™ QmfeO • (2-21) 

Both Eik and Hn- , so defined, are the components of 3-dimensional (spatial) 
tensors, obeying the algebraic restrictions (symmetric, traceless): 

Eik = i^(ife) , h^'' Eik = ; (2.22) 

Hik = -f^(ifc) , h^'' Hik = , (2.23) 

where round brackets denote symmctrisation. By analogy with the vacuum 
Maxwell case for Ei , Bi above, here Eik and Hik also obey differential con- 
straints on the bounding 3-surface. From the Bianchi identities [42,43], one has 
{in vacuo) 

^VkE'^ = , 3y^^ife ^ , (2.24) 

where "^V^ denotes the intrinsic 3-dimensional covariant derivative, which pre- 
serves the 3-metric hij . 

The classical Einstein-Hilbcrt action functional for gravity, with magnetic 
data Hik specified on the boundaries, will be discussed in Sees. 5, 7-9 below, for 
the case of weak anisotropic perturbations. 

3 Boundary conditions in two- component spinor 
language 

As mentioned in the Introduction, a more unified view of the boundary condi- 
tions for perturbed data, as specified on the initial and final space-like hypcrsur- 
faces E/ and E f , can be gained from their description in terms of 2-component 
spinors [15,42,43]. In this Section, we again begin with s = 1 Maxwell pertur- 
bations. 

Consider (in Lorentzian signature) a real Maxwell field-strength tensor F^i, , 
obeying the Maxwell equations (2.2,4). In the theory of 2-component spinors, a 
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space-time index /i is related to a pair of spinor indices AA' { A = ,1 ; A' = 
0' , 1') through the (hermitian) spinor-valued l-forms e'^'^^ , defined by [15,42,43]: 

e^^; = e'^^a,^^' . (3.1) 

Here, denotes a (pseudo-)orthonormal basis of l-forms (a = 0,1,2,3), 
while cr^"^^ denotes the Infeld-van der Waerden translation symbols [15,42,43]. 
At this point, one has to make a definite decision between SL{2,C) spinors, 
which are the most appropriate for studying real Lorentzian-signature geometry, 
and 50(4) spinors, in which Riemannian geometry is most simply described 
[54]. We follow the Lorentzian spinor conventions of [15]. This does not prevent 
one from describing complex or Riemannian geometry - one simply allows the 
'tetrad' of l-forms above to become suitably complex. 
Knowledge of F^^, at a point is equivalent to knowledge of 

Faa'BB' = F^,^ e^^, ^ egg, " (3.2) 

at that point. Here, Faa'BBi is hermitian, for real F^u . Further, the space-time 
antisymmetry F^^ = implies that the decomposition 

FaA'BB' = SAB (j^A'B' + <^A'B' (f'AB (3.3) 

holds, where e^s and ca'B' are the unprimed and primed alternating spinors 
[42,43], and where 

^AB = \Faa'B^' = <I>ba (3.4) 

is a symmetric spinor. In the present Lorentzian case with real Maxwell field, 
(pA'B' is the spinor hermitian-conjugatc to 4)ab ■ (In the Riemannian context, 
Eq.(3.3) would give the splitting of the Maxwell field strength into its self-dual 
and anti-self-dual parts [26,44].) Knowledge of the 3 complex components of 
(pAB at a point is equivalent to knowledge of the 6 real components of Fp,i, at 
that point; also, the <f)AB are, in principle, physically measurable, just as the 
F^v are. 

In terms of the dual field strength *F^v of Eq.(2.5), one finds that 

F^,. + i*F^, = 2<^^Be^'B'e^^>^^', , (3.5) 

together with the conjugate equation. The vacuum Maxwell field equations 
(2.2,4) can then be combined to give 

V^^Vab = , V^^'^A'B' = , (3.6) 

where V^^' = e^^'^^ . 

Here, we are again interested in the decomposition of the Maxwell field 

strength with respect to a space-like bounding hypcrsurfacc and its associated 
unit (future-directed) normal vector n'^ . Define the normal spinor 

n^^' = n^e^^\ . (3.7) 
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Then the (purely spatial) electric and magnetic field vectors and Bk can be 
expressed through 

Ek + iBk = 2<j)ABn'^B'e''''' k , (3-8) 

Eu~iBk = 2 4>A'B' Ub^' e^^'^ . (3.9) 

In 4-vector language, the corresponding co-vector fields E^ and B^ are defined 
by 

= nT^^ , B^ = n''*F,^ , (3.10) 

obeying 

n^'E^ = = nf^B^ . (3.11) 
Next, for e = ± 1 , define 

= 2enVn^B,^^'^'+ . (3.12) 

Here, is symmetric on A and B ; this spinor may be re-expressed in terms 

of Ek and Bk , on making use of the symmetry of rig^ e^^ on its free spinor 
indices B'A' [15]. Here, we define 

.BA'k ^ ftfc^e^^', , (3.13) 



e 



where h'^^ is the inverse spatial metric. The above symmetry property then 

reads 

n^B'^BA'k ^ ns^'e^^''^ . (3.14) 
From Eq.(3.14), we find the decomposition 

-,AB B „AB'k 



{€-l)Ek-i{e + l)Bk . (3.15) 



In particular, our boimdary condition in Sec. 6 below of fixing the magnetic field 
(a spatial co- vector field) on each of the initial and final space-like hypersurfaces 
S/ and Si? is equivalent to fixing the spinorial expression 

= - 2inVe-^^''=Sfc (3.16) 

on each boundary. Note that, even though we regard Bk as having 3 real 
components, the left-hand side, being symmetric on (AB) , appears to have 3 
complex components. In fact, ^'^^ as defined through Eq.(3.16) obeys a further 
hermiticity requirement, appropriate for spinors in 3 Riemannian dimensions 
(that is, on the hypersurfaces E/ and Tip) [15,42,43], so re-balancing matters. 

For comparison with much of the work done on black holes and their pertur- 
bations, one needs the Newman-Penrose formalism [55] - an essentially spinorial 
description of the geometry. Here, considering only unprimed spinors at present, 
a normalised dyad {o^ , i"^) at a point is defined to be a basis for the 2-complcx- 
dimensional vector space of spinors co^ at that point, normalised according to 

OAi^ = 1 = - M . (3.17) 
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The unprimed field strength (f)AB = 4>{AB) can be projected onto the dyad, to 
give the 3 Newman-Penrose quantities 

<?^0 = (I)ABo'^0^ , (1)1 = ^(j)ABO^L^ , (j)2 = (l)AB , (3.18) 

each of which is a complex scalar field (function). Using the Newman- Penrose 
formalism to describe perturbations in the background of a rotating Kerr black- 
hole geometry, Teukolsky [56] derived decoupled separable equations for the 
quantities (j)o {s = 1) and </>2 (s = —1) (for further review, see [39,57].) 
In our non-rotating case, with spherically-symmetric background, the Newman- 
Penrose quantity of most interest to us, following the work of this paper, is 
. In the language of [55,58], has spin and conformal weight zero. This is 
best described in the Kinnersley null tetrad [59] for the Schwarzschild or Kerr 
geometry, in our coordinate system. A null tetrad [55] £^ , , , ifi'^ of 
vectors at a point is a set obeying, in an obvious notation, = n.n = m.m = 
fh.fh = ; i.n = — m.fh = 1 ; £.m = £.fh = n.m = n.fh = . Knowledge of 
such a null tetrad is equivalent to knowledge of the corresponding normalised 
spinor dyad (o^ , t'^) , through the relations 

In terms of the Regge- Wheeler variables to be used for the decomposition of the 
linearised Maxwell field strength F^]) given in Sees. 4, 6, one has 

= i:^T.{^iL + i4t)yemm , (3-20) 

Cm 

where the Yimi^) are the normalised spherical harmonics of [60]. As a result, 
one finds that (/>i obeys the wave equation (4.47) below. With regard to 
the boundary conditions on S/ and Tip , when the variables V^wL '^Um 
being used, the correct boundary data (Sec. 6) will involve specifying both ipuL. 

and ^tV-'ifon on E/ and . 

Turning again to the gravitational field, in 2-component spinor language 
[42,43], one has the decomposition of the Weyl tensor: 

Caf3-fS £AB eCD '^A'B'C'D' + ^A'B' ^C'D' ABCD , (3-21) 

where6 

^ABCD = ^(ABCD) (3.22) 

is the totally symmetric (complex) Weyl spinor, and (in Lorentzian signature) 
^A' B'C D' is its hermitian conjugate. The dual of the Weyl tensor is defined as 

*Cal3'rS = -Z Va/Spa- ^"^^,5 • (3.23) 



(3.19) 
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One finds that: 

(Ca0j5 + i*Ca/3iS^ ^ £A' B' D' ^ ABCD , (3.24) 

together (in Lorentzian signature) with the hermitian-conjugate equation. (If 
instead we had used the Euclidean definition of spinors, then Eq.(3.21) would 
describe the splitting of the Weyl tensor into self-dual and anti-self-dual parts 
[44,58].) The (vacuum) Bianchi identities read [42,43]: 

y^^'^ABCD = , V^^' ^A'B'C'D' = . (3.25) 

On the bounding surface T,f (say), one finds analogously that 

Em +iHu = 2 ^ABCD {n^B' e'^'^'k) {n^ d' e"''' e) (3-26) 
and its hermitian conjugate. Thence, the magnetic tensor Hki is given by 

Hu = (^-i*ABCZ3(n^B'e'''''fc)("Ve^'''^)^ + [h.c.] , (3.27) 

with a corresponding equation for E^^ . These two equations can straightfor- 
wardly be inverted to give an expression analogous to "i^^ for the s = 1 Maxwell 
case of Eq.(3.16). This analogous expression, i^^BCD ^ ^,{^bcd) ^ again to- 
tally symmetric on its indices, and is given for e = ± 1 by the generalisation of 
Eq.(3.12), as 

^ABCD ^ 4 e A, n'' n^' c n"" D' ■ (3-28) 

Thus, again for the magnetic case e = + 1 , ^/^bcd pj-Qvidcs a spinorial version 
of the magnetic part Hk£ of the Weyl tensor, to be fixed on the final bound- 
ary Si? in our s = 2 treatment of the gravitational boundary-value problem, 
perturbed about spherically-symmetric collapse. Of course, the pcrturbative 
boundary data Hki must further be chosen such that the divergence condition 
(2.24) holds: ^VkH^'' = on Sf , just as, from Eq.(2.13) in the Maxwell case, 
the condition ^VkB^ = must hold. 

Two of the five complex components of ^ abcd are contained in the Newman- 
Penrose quantities [42,43,55] 



*o = -^abcdo^o^ 0° , ^-4 = ^abcd i^i'^ i° , 

(3.29) 

where (o^ , t^) is a normalised spinor dyad [43,43,55]. Again taking the Kinner- 
sley null tetrad [39,59], it was further shown by Teukolsky [56] that ^'o and ^'4 
each obey decoupled separable wave equations. Following work by Chrzanowski 
[61], it was confirmed by Wald [62] that, given a solution "^/q or ^'4 of the 
Teukolsky equation for a (nearly-) Kerr background, all the vacuum metric per- 
turbations can be reconstructed in a certain gauge through a series of simple 
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linear operations on ^'o (or on ^'4) [39,61,62]. Once the linearised metric per- 
turbations are known, then, of course, one can compute other Newman- Penrose 
quantities at linearised order, such as 

*2 = ^ABCDO^O^ L^" . (3.30) 

Note here that, for an unperturbed Schwarzschild background (say), only the 
middle Newman-Penrose quantity , out of the set (i = , . . . , 4 ) , is non- 
zero, with ^2 = — M/r^ in Schwarzschild coordinates [57]. 

The analogous process is implicit in the (Maxwell) discussion above, leading 
to Eq.(3.20): For s = 1 perturbations of the Kerr metric, the correspond- 
ing linearised Maxwell vector potential (in a particular gauge) can be re- 
constructed by simple steps from the Newman-Penrosc quantites ^0 or (f>2 
[39,61,62]. Hence, the middle Newman-Penrosc quantity (pi can also be found, 
leading to Eq.(3.20). 

In the more complicated s = 2 case, although we have not yet finished de- 
tailed calculations on this point, it does now look reasonable to expect that, 
for gravitational perturbations about a spherically-symmetric gravitational col- 
lapse, there should exist a relation analogous to the s = 1 Eq.(3.20). In this 
case, for a Schwarzschild background, this relation would involve expanding out 
the first-order perturbations of the middle Newman-Penrose quantity ^'2 in 
terms of the s = 2 Regge- Wheeler description of Sees. 5, 7-9 below. 

4 Regge- Wheeler formalism — Maxwell case 

A more unified treatment of the angular harmonics appearing in the separation 
process for s = 1 (Maxwell) and s = 2 (graviton) perturbations of a spherically- 
symmetric background can be given in terms of vector and tensor harmonics 
[63]. In [11], we expanded the s = (massless-scalar) perturbations in terms 
of scalar spherical harmonics Y£m{9 ,4') ; which have even parity. Vector and 
tensor spherical harmonics, however, can have odd as well as even parity. 

Any vector field in a spherically-symmetric background, such as the classical 
s = 1 (photon) solutions appearing in the background gravitational-collapse 
geometry , can be expanded in terms of vector spherical harmonics on the unit 
2-sphere [37,38]. Correspondingly, angular vector and tensor indices are raised 
and lowered with the metric %b , given by 

7ee = 1 , %<!> = sin^^* , 7(?0 = %e = . (4.1) 

The even-parity vector harmonics [37,38] have angular components 

{^em)a = (daYem) , (4.2) 

where a = {6 ,(j)) ■ The odd-parity vector harmonics are 

($€m)a = ea'(56^€m) • (4.3) 
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Here, e^j denotes the tensor with respect to angular indices (a = S ; b 
9 ,(f)) , such that the lowered version Cab = — eta is anti-symmetric, with eoi 

(7) ^ =syn9 , where 7 = det(7o6) . Thus, 



-1 



= smt 



= 



(4.4) 



The forms of the angular harmonics appearing in the s = 1 photon calculations 
below (;:an be deduced from these vector-spherical-harmonic expressions. 

Analogously, any rank-2 tensor field such as a linearised (graviton) s = 2 
classical solution, to be treated in Sees. 5, 7-9, below, can be expanded in terms 
of tensor spherical harmonics. The even-parity tensor harmonics are 



{^tm)ab = Yi, 



,\ah 



{^im)ah = Jab V^ri 



(4.5) 



where a bar | denotes a covariant derivative with respect to the metric ■jab ■ The 
odd-parity tensor harmonics are 



(X^m)ab 



(4.6) 



Prom 1957, Regge and Wheeler [63] developed the formalism for treating 

both spin-1 and spin-2 classical perturbations of the Schwarzschild solution and 
of other spherically-symmetric solutions, corresponding to Maxwell and gravi- 
tational (graviton) perturbations. Here, in the Regge- Wheeler (RW) formalism, 
for Maxwell theory we decompose the real linearised field strength pji]) and lin- 
earised vector potential A^j^ '' into tensor and vector spherical harmonics, respec- 
tively [37,38] . We are assuming that the background (unperturbed) classical so- 
lution consists, as above, of a spherically-symmetric gravitational and massless- 
scalar field (7^,^ , $), with no background Maxwell field: aI^*' — , F^"' — . 

For each spin s = , 1 , 2 , the corresponding perturbation modes split into 
those with even parity and those with odd parity. Under the parity inversion: 
^ — » {tt — 9) , (j) ^ {tt + (p) , we define the even perturbations as those with 
parity tt = (—1)^ , while the odd perturbations have parity tt — (— l)^"*"^ . For 
Maxwell theory (s = 1) , the £ = mode corresponds to a static perturbation, in 
which a small amount of electric charge is added to the black hole; in particular, 
a Schwarzschild solution will be 'displaced' infinitesimally along the family of 
Reissner-Nordstrom solutions. For radiative modes with £ = 1 (dipole) and 
higher, we set 



00 e 



F'^Xx) = EE (^i^ 
= E E 



=1 m=-e 



em 



em 



em 



em 



(4.7) 
(4.8) 
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On substituting this decomposition into the boundary expression (2.7) for the 
classical Maxwell action S^i^^ , we find 



t^m'J Jo V Jem\ Jem' 



(4.9) 

where Q. denotes the angular coordinates ^ , </> , and we work at present with 
the Lorentzian form of the spherically-symmetric background metric (for which 
the Riemannian form was given in Eq.(1.7), by suitable choice of coordinates), 
writing: 

= - e''^*''-) dt^ + e«(*''') dr^ + {dO'^ + sin^ e dcj)"^) . (4.10) 

For later reference, we introduce a standard expression for such a spherically- 
symmetric space-time, defining the 'mass function' m{t , r) by 

e-« = 1 - , (4.11) 

r 

within the Vaidya-likc region of the space-time, in which the blaek hole is evap- 
orating. Clearly, the odd and even contributions decouple in the action of 
Eq.(4.9). 

For the subsequent detailed treatment of the angular harmonics involved, 
we follow Zerilli's decomposition [64] of and ^^^^ . Thus, we set 

em 

I h /2 h \ 

V-/3 -h -h / 

For a given choice of (^ , m) , we take 

h = (rt^) Ye^m , (4.14) 
f, = (sm^?)-V/^(;)) {d^Y,„,) + (fI^A (a«y,„0 , (4.15) 

\ / im \ / 1771 

\ / Im \ / tm 

U = {smOyUF^^A (d^Yern) + (F^!:f) {deYe,n) , (4.17) 

\ / im \ / im 

h = - (Flf) {&'^n6){dgY,ra) + (Plf) {d^Ytm) , (4.18) 
\ / im \ J im, 

= (F(J)^^(sin0) y,„(17) . (4.19) 
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Here, the Yfm{^) are scalar spherical harmonics [60], and a caret indicates that 
the quantity is a function of t and r only. 

Again, following [63], for the vector potential we set 



/^(o)\ _ Q a2em{t,r){d4,Yim) - a2im{t,r) (sine) (deYgrn) 



(^+I)sin6» ' £{i+l) 
- aoemit, r) Yfm(O) , aumit, r) Yim{Q) ,0,0 



4.20) 
(4.21) 



Now Eq.(4.9) can be expanded out in the form: 



1 n pRoO 



-6)/2 



e'm' 
Roo 



(4.22) 



Of course, the components of the field strength are given in terms of those 
of the vector potential by Eq.(2.3); for example, F^l^ = dtA^^'' — deA^^^ . This 
gives the relations 



t9 

pie 
tr 

{o 



Km 

im 

Im 
Im 



{dta2lm) 

Holm , 
{dra2im) 

{drdoem) - (dtaum) 

0'2em ■ 



(4.23) 
(4.24) 
(4.25) 
(4.26) 
(4.27) 
(4.28) 



The action (4.22) then simplifies to give 



^class 



1 i-noo 



where 



f(o) _ 
Jim ~ 



i{i + 1) 



(4.29) 

(4.30) 
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determines the odd-parity Maxwell tensor via Eqs. (4.23,25,28). 

The form of the classical action (4.29) can be further simplified by using the 
Maxwell field equations (2.2,4). This will lead finally to the form (4.45) below, 
in which 5*^^^ is expressed explicitly in terms of bomidary data, as needed in 
the subsequent calculation of the quantum amplitude (see [11,16] for the spin-0 
analogue). 

The linearised Maxwell equations can be written as [49] : 



The fi 



t , r equations give 



= 



dr 



^(£+l)e»(F(;)) 



Im 



im 



dr 



e-"£(£ + l) 



im 



\ / im 



(4.31) 

(4.32) 
(4.33) 
(4.34) 
(4.35) 



The ^ = , 4> components give the same equations. Note that Eq.(4.32) is just 
the (source-free) constraint equation di£^^^^ = of Eq.(2.14). 

The equations (4.23,25,28) together imply the decoupled wave equation for 
odd perturbations: 



where 



Vu{r) = 



Vu{r) a2i 

£+1) 







> 



(4.36) 



(4.37) 



is the (massless) spin-1 effective potential and where, as usual [63], we write 
dr' = e~° dr ■ If the geometry were exactly Schwarzschild, then the coordinate 
r* so defined would be the Reggc- Wheeler or 'tortoise' radial coordinate, given 
by [49,63]: 

r* = r + 2Mln(r-2M) . (4.38) 

As in [11,16], we assume that the adiabatic approximation is valid in a neigh- 
bourhood of the initial and final surfaces, S/ and • In that case, we can, as 
before, effectively work with the field equations on a Schwarzschild background, 
except that the Schwarzschild mass Mq is replaced by the mass function m{t ,r) , 
as defined in Eq.(4.11), where m{t , r) varies extremely slowly with respect both 
to time and to radius. 

Equation (4.33) gives dt{F'^g^)im also in terms of {F^''g ))>m , while Eq.(4.35) 



gives dt{F^r^)ira also in terms of {F^f) 
that 



Im ■ Together, Eqs. (4.27,33,35) imply 
= , (4.39) 
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where we define 

&2 = e-'^aum . (4.40) 

Thus, with a suitably defined variable , the even perturbations obey the 
same decoupled wave equation (4.36) as the odd perturbations. 
Finally [65], we set 

2 

4mit,r) = j^f^^({dtaum) - (draoem)) , (4.41) 

which is clearly gauge- invariant. Now V'^^ is related simply to the (even-parity) 
function of the previous paragraph: Eqs. (4.33,39) imply that 

{dt4':l) = - /iri • (4.42) 

Equations (4.40,41) can now be used to simplify the classical Lorcntzian 
action (4.29). For ease of comparison with the (second- variation) classical spin- 
2 action, where the pattern is similar, we define 

4tn = «2^m(*.0 , (4-43) 

4t = i{i+l)4tit^r) . (4.44) 

Then, given weak-field Maxwell boundary data, specified by the linearised magnetic- 
field mode components on each of the boundaries S/ and Si? , the 
corresponding classical Maxwell action is 

(4.45) 

Of course, the limit i?oo oo must be understood in Eq.(4.45). 
Note further that, from Eqs. (4.36,39), one has 

(a„ = - aaem , (4.46) 

whence ^^''^ obeys the same decoupled wave equation (4.36,39) as for 

a2em and for /j^ , namely, 

i9r*f4t-i9tr4t-Vu4t = • (4-47) 

The spin-1 radial equation (4.36,39,47) in a Schwarzschild background, both for 
odd- and even-parity Maxwell fields, was first given in 1962 by Wheeler [66]. 

This suggests a 'preferred route' for understanding the even-parity pertur- 
bations (which are more complicated than in the odd-parity case, which only 
involves the single function a2em{t,r) , obeying the decoupled field equation 
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(4.36)): Given suitable boundary conditions, one first solves the linear decou- 
pled wave equation in two variables t and r , namely, Eq.(4.47), for ipumi'^ y ^) ■ 
By differentiation, following Eqs. (4.44,46), one then finds aoem{t , r) . Then, by 
integrating Eq.(4.41), one obtains also ai^„i(t ,r) , and hence, from Eq.(4.40), 
ftm ' ■ Fi'oni Eqs. (4.23-28), one now has all the non-zero components of the 
Maxwell field strength F^i, in this even-parity case. 

5 Regge- Wheeler formalism — odd-parity gravi- 
tational perturbations 

Our boundary- value problem, as posed in the Introduction and in [11,16], in- 
volves specifying on the final space- like hypersurface Y^p the spatial compo- 
nents h^j^\x) {i,j = 1,2,3), of the real perturbations h^i}}{x) of the 4-mctric 
= 0,1,2,3). We shall construct a basis of tensor spherical harmonics 
with which to expand the angular dependence of /i-j^ . In general, we make a 
multipole decomposition for real metric perturbations, of the form: 

oc i 

f^i^) = E E + ' (5-1) 

£=2 m=-e 

where — and -|- denote odd- and even-parity contributions, respectively. We 
comment below on the limit = 2 in the summation over £ . 

In Sec. 4, s = 1 (Maxwell) perturbations of spherically-symmetric gravita- 
tional backgrounds were treated in the Regge- Wheeler (RW) formalism [63], 
and split naturally into odd and even type, according to their behaviour under 
parity inversion: 9^ {tt — 9) , (f) ^ {-k + (()). The even 'electric-type' perturba- 
tions have parity tt = (—1)^ , while the odd 'magnetic-type' perturbations have 
parity tt = (—1)^+^ . The analogous (orthogonal) decomposition also holds for 
the s = 2 gravitational-wave perturbations. 

In the s = 1 Maxwell case, the lowest t — Q mode does not propagate: 
in the electric case, it corresponds to the addition of a small charge to the 
black hole, to turn a Schwarzschild solution into a Reissner-Nordstrom solution 
[49] with charge Q<^M . Correspondingly, in the s = 2 case of gravitational 
perturbations, the multipoles with £ < 2 are non-radiatable. For example, 
the even-parity gravitational perturbations with C = correspond to a small 
static charge in the mass, while the £ = odd-parity perturbation is identically 
zero. For £ = 1 , the odd-parity (dipole) gravitational perturbations must be 
stationary [39] , and even-parity dipole perturbations correspond to a coordinate 
displacement of the origin [67] and can be removed by a gauge transformation. 
For a general spin s = 0,1,2, perturbations with £ < \s\ relate to total 
conserved quantities in the system. In the present s = 2 gravitational-wave 
case, we consider accordingly only the propagating £ — 2 (quadrupole) and 
higher-.^ modes. 
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In this Section and in Sec. 7 below, we restrict attention to odd-parity gravitational- 
wave perturbations. Following Moncricf [68], we write 



(5.2) 



(N.B. : one should not confuse the subscripts 1,2 here with spin subscripts.) 
The non-zero components of the symmetric tensor fields [(ei 2) jj] ^ are defined 
by 

[(ei)re]^^ = - {d^Yi^) / [sme) , 
[(ei)r<^]^„ = {sine){d0Yem) , 



1 

2 



[(62) 

[{e2),j,4>)]im 

This basis is normalised according to: 



= (sin e) ^ [(sin 9) d^^ - (cos 0) d, 
(sin 6*)"^ dl - (cos 6*) de - {sm.6) 



Yen 



(5.3) 
(5.4) 

(5.5) 



(cos^) de - {sm6)dg^ 



Yen 



, 
26-" 



i{l+l)5ee'5mrr. 

2r4 



5ee' S„ 



(5.6) 
(5.7) 

(5.8) 
(5.9) 
(5.10) 



where /dfi ( ) = j'^^ d<l) d6 smd ( ) , and where these indices are raised 
and lowered using the background 3-metric 7*-' ,7ij . Note that both [(ei)^]^^ 
and [(e2)ij]^^ are traceless. 

In the standard 3 + 1 decomposition for the gravitational field [49], the 4- 
metric g^^, is decomposed into the spatial 3-metric hij = gij on a hypersurface 
{x^ = const.} , together with the lapse function N and the shift vector field , 
such that the 4-dimensional space-time metric has the form [15,49]: 



ds'^ = hij {dx' + N' dt) {dx^ + dt) - N'^di^ 
For odd-parity perturbations of the lapse, one has 

7v(i)(-) = , 

while the odd-parity shift vector takes the form 



(-) 



em 



0, 



1 



(sin6») 



),{sme){deYem) 



For a real 4-metric g^i, , both h[^^ and N^^''^ are real, and one has 
/,(-)* _ (_-[\m l(-) 



(5.11) 
(5.12) 

(5.13) 
(5.14) 
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In the Hamiltonian formulation of general relativity, the momentum tt*-' 
conjugate to the 'coordinate' hij is a symmetric spatial tensor density [49]. As 
with the 3-metric hij above [Eq.(5.1)], the linearised perturbations of Wij can 
be decomposed into multipoles with odd or even parity: 

oo e 

^lj'(^) = 1212[i^ll')eJ^) + i^il')eJ^)] ■ (5-15) 
e=2 m=-e 

Restricting attention at present to the odd modes, one has 

= i^'h)^ {Pumit,r)[{e,),,]^^+p2e^{t,r)[{e2h],^) , (5.16) 
where [(ei)jj]^^ and [(e2)ij]^^ are given above. One finds that 

Pum{t,r) = ^ ((dth[-l] -r^dri'-^]] , (5.17) 




2N(o) 

P2im{t,r) = -^^l^(d,h^-l)+2h[-l}l . (5.18) 

One can typically simplify the form of the perturbations by a gauge transfor- 
mation (linearised coordinate transformation) in a neighbourhood of the final 
space-like hypersurface S p . Suppose that the infinitesimal transformation is 
along a vector field . Then the metric perturbations transform infinitesi- 
mally by 

Qiiv g,i,y - VpCi^ - V,y£,f, . (5.19) 

For odd perturbations, in the notation of Eq.(4.2), consider the infinitesimal 
vector field , with components [63] given by 







(5.20) 



The resulting 'gauge transformation' is summarised by 

= hi-^l-{dtAim) , (5.21) 

h[-Z = h[-l-{drAim) + ^ , (5.22) 

h^2~Z = h[~l + 2K,^ . (5.23) 

We have here neglected time-derivatives of the metric components: we are 
assuming that an Ansatz for the gauge functions A£„i(t , r) based on separation 
of variables will be valid, involving frequencies which satisfy the adiabatic ap- 
proximation [11,16] described below. In the Regge- Wheeler gauge [63], we set 

^olrif^ = ^otli ^^'^ = ^uli ' ^ Eqs. (5.21,22), but require also 

h'k^'^ = = + 2 A,™ . (5.24) 
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For each £ , one can obtain solutions for arbitrary m by rotation from the case 
m = . Note also that the above equations show how the RW perturbations 
can be uniquely recovered from the perturbations in an arbitrary gauge. 

Since odd-and even-parity perturbations decouple, the odd-parity field equa- 
tions in the RW gauge are obtained by substituting Eq.(2.7), together with the 
equation N-~^ = h[^^ for the linearised shift vector, into the source-free lin- 
earised Einstein field equations [49]. The (Lorentzian) spherically-symmetric 
background metric is taken, as above, in the form of Eq.(4.10), and the mass 
function m{t,r) is defined by Eq.(4.11). Then the odd-parity linearised field 
equations, taking respectively the {t(j)) , {rcj)) and {9(1)) components, read: 



" - dAh^uT - 9^h[-,r + Fu{t , r) h[^r = , (5.25) 

and 



+ K« + ^)(^'-^-if"-5*/.l;r-) , (5.26) 

i(a + 6)+i6'(a-6))47if--F..(i,r)4if- = 

5*(e^"-^'/^472.^^)-9.(e(^-<^)/^/^i72.^^) = . (5.27) 

Fu(t,r) = ^ + 4m' - £(^ + 1)^ + e"-''(^a + i o(d - fe)^ + Z , (5.28) 



and 
Here 



and 



pb 1 

F2i{t,r) = - ^{£ + 2){£-l) + Ze^-''+ d+-a{a-b) , (5.29) 



with 

(m' + rrnj | . (5.30) 



2e" / ,, f2m'e'' 
Z = \m"+( — ^ 

Further, the Einstein field equations imply 

m' = Anr^p . (5.31) 

Here, p is the total energy density of all the radiative fields; in the present case 
p has contributions from s = (massless scalar), s = 2 (graviton) and, if the 
Lagrangian contains Maxwell or Yang-Mills fields, also s = 1 . In the adiabatic 
limit, in which m{t,r) varies extremely slowly, Eq.(5.31) provides part of the 
description of the approximate Vaidya metric for the region of space-time in 
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which the black hole is evaporating; a much fuller treatment of the Vaidya 
region is given in [69]. 

Now consider the lowest-order perturbative contribution to the classical 
action of our present coupled bosonic system, for which the non-zero back- 
ground part consists of the spherically-symmetric gravitational and scalar fields, 
(7/iL/ 5 ^) • It is assumed that the background fields obey the coupled Einstein/massless- 
scalar field equations and contribute S'^J'^gg to the classical action. Any spin-1 
fields present, whether Maxwell or Yang-Mills, propagate at lowest order in the 
background, and the spin-1 contribution to the classical action begins at second 
order, as described in Sec. 4, and adds to the various lowest-order perturbative 
Einstein-Hilbcrt and scalar contributions. To simplify the exposition at this 
point, let us temporarily neglect the spin-1 field. 

After detailed calculation [33,34], one finds that the classical Lorentzian ac- 
tion for an (anisotropic) Einstein/massless-scalar solution, subject to perturbed 

boundary data (h^j^ + /i-j"* , $ -|- f/)^^^) on the initial and final surfaces S/ 

and S F , can be written as the background contribution S^J^^^^ above, plus a 
quadratic-order contribution, plus higher-order terms. Thus: 



Again, for simplicity, it is assumed that no perturbations arc prescribed on the 
initial boundary S/ , but that there are non-zero perturbations on the final 
boundary T,f ■ As described after Eq.(5.14), tt'-' is defined to be the momentum 
conjugate to the 'coordinate' hij . Similarly, tt^ is the momentum conjugate to 
the variable 4> ) and is given by the normal future- directed derivative d4>/dn . In 
Eq.(5.32), T, as usual, denotes the Lorentzian proper-time interval between the 
initial and final boundaries, as measured at spatial infinity. Further, M denotes 
the mass of the space-time; for a well-posed (complexified) classical boundary- 
value problem, the mass defined on the initial surface must agree with the mass 
defined on the final surface. 

We can now compute the odd-parity contribution to the classical gravita- 
tional action. In an arbitrary gauge, taking only the intrinsic metric perturba- 
tion h^^^ to be non-zero on Hp , but (f)^^'> = there, one has the second-variation 
part of the action 



On discarding a total divergence, the spin-2 classical action can also be written 



o, _ , 0(2; 

•Jclass — '-'class ^ '-'cla 




- MT 



(5.32) 




(5.33) 
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as 




(5.34) 



Since odd- and even-parity perturbations are orthogonal, there are no cross- 
terms in the action. In an arbitrary gauge, the odd-parity contribution to the 
classical gravitational action can be re-written as 




(5.35) 

Here, wc have used Eqs. (5.8-10,17,18), and have taken the perturbations to van- 
ish initially. Note that, if we were to evaluate Eq.(5.35) in the RW gauge, for 
which h2il^^ = , so that the second integral would vanish, and then substi- 
tute for the Regge- Wheeler functions via Eqs. (5.21-24), then, in the adiabatic 
approximation, we would arrive back at Eq.(5.35) up to a boundary term 

hi'eL* PL^r'''^ , (5-36) 

where 

Thus, S^fl,^[{h\j^)im] is gauge-invariant up to a boundary term. We shall 
return below to the question of boundary conditions for the odd-parity pertur- 
bations. 

At first sight, the odd-parity action looks unwieldy. Ideally, we would like 
to work with a classical action (both for odd and for even parity) in the form 
/ drip{dtil^) , of the same general kind as in the massless-scalar classical action 
in [33,34]. In the present gravitational case, V would ideally also be gauge- 
invariant and would obey a wave equation with a real potential. To achieve this 
form, first use Eqs. (5.21-24) for the RW functions, and substitute them into the 
field equation (5.26), to obtain 

{dtPL^) = n£+miir)U-l + ldr(f^yj , (5.38) 
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{drPL^) = - ^ + ^(^ + 1) (iFuir) + ^) ({dthi-l) + k39) 



When we substitute into Eq.(5.35), using Eq.(5.38) for h[^^l^, and then use 
Eqs.(5.29,39), the boundary term (5.36) vanishes. As a consequence, we find in 
the adiabatic approximation that 



o(2) 
class 



where ^2im defined by 



(5.40) 



41. = rP^'' ■ (5.41) 



Eqs. (5.21-23) show that ^2im is gauge-invariant. Indeed, ^2em is related 



to Moncrief's [68] gauge-invariant generafisation of the ZerilH function, ^ 
[64,70], defined as 



by 

Q^fJ = - 1^ (5*41) , (5.43) 

which, in effect, replicates Eq.(5.38). Note the simplifying property of Q^"^ , 
namely, that it is written entirely in terms of perturbations of the 3-geometry 
(our chosen boundary data). Further, <3^~^ is automatically gauge- invariant, 
since it is independent of the perturbed lapse and shift. 

In the adiabatic approximation, the function obeys the wave equation, 
of RW type: 

e-«9.(e-«(a.41)) - {dtf^ijl - vt\r)Ct^ = , (5.44) 



where 



Further, in the adiabatic approximation, the function Q'"^^^ obeys the same 
equation (5.44). In the RW gauge, one would solve for Q[^^^ , then determine 
^Urn^^ from Eq.(5.42), and then determine h^^^^ with the help of Eq. (5.27). 



6 Boundary conditions and classical action — Maxwell 
case 

In Sec. 4, for the Maxwell field (whether odd- or even-parity), we derived an 
expression, Eq.(4.45), for the classical Maxwell action Sf^^ , given via certain 
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operations in terms of the magnetic field on the final boundary T,p . In Sec. 5, 
for odd-parity gravitational perturbations, we derived Eq.(5.40), a somewhat 
analogous expression for the second- variation classical action S^^^^^ , depending 
on the odd-parity metric perturbations h^J^ over the final surface Si? . In the 
present Section 6, concerning Maxwell theory, we relate Sf^^ more explicitly 

(2) 

to the final boundary data. The following Sec. 7 will similarly re-express S^^^^^ 
for odd gravitational perturbations. A corresponding treatment for the more 
complicated even-parity gravitational case will be given in Sees. 8, 9. The result- 
ing expressions for the classical action will greatly simplify our understanding 
of the black-hole quantum amplitudes. 

Physically, our gauge-invariant odd- and even-parity Maxwell variables > 
V'wL ^I'c effectively the radial components of the magnetic and electric field 
strengths, respectively: 

= ^4tit,r)Yerr.m , (6.1) 
^l^'^(^) = - ^V'[l(i,r)r,„(17) , (6.2) 



where 

B^H' = (^'V)"' • (6.3) 

The remaining, transverse, magnetic field components are 



Blll'ix) = 



-a/2 



-a/2 



(sin 0) 



aumV',r)[d0Yim) H 



'+1) (sir 



(6.4) 
. (6.5) 



The main aim of this paper is to calculate quantum amplitudes for weak 
spin-1 (Maxwell) perturbative data on the late-time final surface Si? , by eval- 
uating the classical action S'^jjfg and hence the semi-classical wave fimction 
(const.) X exp(i 

■^ciass)' functionals of the spin-1 final boundary data. In 
Eq.(4.45), S^^g was expressed as an integral over the boundary, involving var- 
ious perturbative quantities used in the description above of the dynamical 
perturbations. The present task is to determine 'optimal' or 'natural' bound- 
ary data, both for the odd-parity case and separately for the even-parity case, 
such that (i.) the classical boundary-value problem can readily be solved, given 
these data, and (ii.) the classical Maxwell action S^^^ can be (re-) expressed in 
terms of the appropriate boundary data. Under those conditions, we will then 
have a description of the spin-1 radiation, associated with gravitational collapse 
to a black hole, analogous to that for the spin-0 (massless-scalar) radiation, 
developed in [33,34]. 

Following the discussion of Sees. 2 and 4, the relevant field components to 
be fixed on S/ and Sj?, in the odd-parity case, are tpuin' ^® ^^^^ 
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from Eq.(6.1). In the even-parity case, the appropriate boundary conditions 
involve fixing aum , as may be seen from Eqs. (6.4,5). In the even case, this is 
equivalent, from Eqs. (4.40,42), to specifying dtipuln boundary. Thus, 

we choose the boundary data to consist of tpuin odd-parity case, and 

^ffpieln ^^'^ even-parity case. Hence, though V'hL and ijj'fdn both obey 
the same dynamical field equations (4.36,43) and (4.47), the natural boundary 
conditions for them on S/^f are quite different - Dirichlet for the odd-parity 

case V'wL ' Neumann in the even-parity case V'hL • "This is reminiscent 
of the situation obtaining when spin-2 gravity is coupled to all lower spins, 
fermionic as well as bosonic; that is, to spins s = | , 1 , ^ and , especially in 
locally-supersymmetric models [22] , such as models of AT = 1 supergravity with 
gauged supermatter [71-73]. As found also, in this paper, for spin-2 (graviton) 
perturbed data on the final surface S f , the natural boundary conditions are 
contrasting, for odd-parity vis — d — vis even-parity modes. For the remaining 
s = (scalar) bosonic case, if one requires, as in the Introduction, that the full 
theory be locally supersymmetric, so that quantum amplitudes are meaningful, 
one finds that all scalar fields must be complex . whether as part of a multiplet 
or as a single complex scalar [22] . The treatment of the real s = case in [33,34] 
can be replicated in the case of a complex scalar field (j) , except that the natural 
boundary conditions, consistent with local supcrsymmctry, require Re((/)) to 
be fixed at a surface such as T,p (Dirichlet), whereas the normal derivative 
9[Im((/))]/9n must also be fixed at a bounding surface (Neumann). Of course, 
this treatment extends to fermionic data (s = ^ and |) , as described in [35,36]. 

In the gravitational-collapse model, by analogy with the simplifying choice 
4>^^^\-^j = for the initial perturbative scalar-field data, taken in [33,34], we take 
(for the purposes of exposition) the simplest Maxwell initial data at S/ (t = 0) . 
That is, we consider a negligibly weak magnetic field outside the 'star': the 
boundary conditions are 



Condition (6.6) is a Dirichlet condition on the initial odd-parity magnetic field 
- see Eqs. (6. 1,4, 5). Condition (6.7) implies that we have an initially static 
even-parity multipolc [65]. 

We now follow the analysis of the spin-0 field, and separate the radial- 
and time-dependence. In neighbourhoods of S/ and Sjr, where an adiabatic 
approximation is valid, we can 'Fourier-expand' the variables V'wm(*''') 
V'hLC^'^) ' subject to the initial conditions (6.6) and (6.7). By analogy with 
the scalar case [33,34], let us write 



V'il(0,r) = 
9*V'll(0,r) = 



(6.6) 
(6.7) 




(6.8) 



(6.9) 
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where the radial functions {ijj[°kf {r)} and {ijj['kf {''')} arc independent of m , given 
the spherical symmetry of the background space-time. Here, the position- 
independent quantities {a^'^^^} and {a['^k£m} some coefBcients, with smooth 
dependence on the continuous variable k , which label the configuration of the 
electromagnetic field on the final surface T, p . 

The radial functions {tpikei''')} and {'ipikei''^)} ^&ch obey a regularity condi- 
tion at the centre of spherical symmetry {r = 0} on the final surface S p ; this 
requires that the corresponding (spatial) electric or magnetic field, defined via 
Eqs.(6.1-5), should be smooth in a neighbourhood of {r = 0} . As a consequence, 
the radial functions must be real: 

= 4%,ir) , = 4%,ir) ■ (6.10) 

For small r, the radial functions should be asymptotically proportional to a 
spherical Bessel function [74]: 

^ifek'') ~ (const.) X rji{kr) , (6.11) 
■^ikei'^) ~ (const.)' X rji{kr) , (6.12) 

as r — > 0+ . Also, the reality of the radial electric and magnetic fields implies 
that 

4tit,r) = i-ir 4t-mit,r) , 4fmit,r) = (-1)" ^t-mit^r) . 

(6.13) 

This in turn implies that 

«S™ = (-I)'" , = (-1)'"+' 4:^;.,-™ • (6.14) 

Since the potential (4.37), appearing in the {t,r) wave equation (4.36), tends 
sufficiently rapidly to zero as r ^ cxd , in the region where the space-time is 
almost Schwarzschild, one has asymptotic (large-r) behaviour of V'ifcl(^) a^^d 
^^[^{r) which is analogous to that in the scalar case: 

4kiir) ~ (4°)exp(ifcr:)+zif exp(-zfcr:)) . (6.15) 
4keir) - [zif eM^kr:) + zit eM-^kr:)) . (6.16) 

Here, {z^^))} and {^fc^''} are complex coefficients, depending smoothly on the 
continuous variable k . Also, r* again denotes the Regge- Wheeler 'tortoise' 
coordinate of Eq. (4.38) for the Schwarzschild geometry [49,63]. As in the scalar 
case [33,34], the inner product (normalisation) for the radial functions follows 
in the limit i?oo — > oo : 

r°°dre«^SWV'lfcUO = 2n\ziff\s{k,k') + 5{k-k')\, (6.17) 

Jo '- ^ 



rRc 

Jo 



dre^4^,]eir)4teir) = 2n\zil^\' \s{k ,k') + 6{k ,-k')] . (6.18) 
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Finally, wc arc in a position to compute the classical Maxwell action S^^^^ 
as a functional of the spin-1 boundary data on the final surface Ei? . This gives 
straightforwardly the semi-classical wave function for the complexified time- 
interval T , leading to the Lorcntzian quantum amplitude or wave function. Our 
boundary conditions (6.6,7) above on the initial hypersurface E/ , at time t = 
say, were designed so as to give zero contribution from S/ to the expression 
(4.45) for the classical action S^^^. The contribution to (4.45) from Si? is 
found, using Eqs.(6.8,9,17,18), to be 

OEM U Jo) (e) 1-1 

'-'class Ll"lfe£m ' "lfc£mJj 

~ 2^(^ + l)!7o [' '^'^ ' y ^keml +^^\"-lkem"'l-kemJ J (g j^g^ 



As promised, this does now express the classical Maxwell part S^^^ of the ac- 
tion as an explicit functional of suitably chosen boundary data, namely, {difc^^} 
and {auLl • 



7 Boundary conditions and asymptotically-flat 
gauge — odd-parity gravitational perturbations 

In classical Lorcntzian general relativity, one would expect to choose regular 
Cauchy data on an initial space-like hypersurface S/ , which would then evolve 
smoothly into {x^ > 0}, subject to the linear hyperbolic equation (5.44). A 
natural initial condition, for given quantum numbers {£,m) [65], would be to 
assume an initially stationary odd-parity multipole; 

= . (7.1) 

The combined Einstein/massless-scalar boundary- value problem, originally 
posed in [33,34], for complex time-separation T = \T\ exp{—i9) , < 9 < Tr/2 , 
involved specifying the intrinsic 3-metric {hij)i^F and the value of the scalar field 
{<P)i,F on the initial and final space-like hypersurfaces S/ , Si? . By Eq.(5.43), 
the above Eq.(7.1) reads 

QL'(0.0 = (7.2) 

or, equivalently, 

h['il{0,r) = , (7.3) 
hi'eliO.r) = , (7.4) 

We therefore take these as our (Dirichlet) boundary conditions on the odd- 
parity gravitational perturbations, on the initial surface S/ , even though they 
may have arisen from consideration of the Cauchy problem. 
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In [11,16,33,34] for the s = case, we made use of the adiabatic approxima- 
tion in order to separate the perturbation problem with respect to the variables 
t and r . Here, for s = 2 , we first separate the odd-parity Eqs. (5.25,26) in the 
RW gauge, and then use Eqs. (5.21-23) to determine the time-dependence (in 
particular) in any gauge. 

As in the massless-scalar (s = 0) case, we introduce a 'Fourier-type' expan- 
sion: 

/oo 
-oo 

where the {dkem} certain odd-parity 'Fourier 'coefficients. With suitable 

treatment of any arbitrary phase factors involved, in order to separate the odd- 
parity field equations (5.22,23) in the adiabatic approximation, one must have 

hMm'^it^r) « cos(fct) , (7.6) 
oc sm{kt) , (7.7) 

(Of course, a normal-mode e~''^* time dependence for the functions h^'l^^ 

would also satisfy the field equations.) In Eq.(5.23), if h'^^^ , which is related 
to the odd-parity shift and can thus be freely specified, has cos{kt) time depen- 
dence, then /io7m must have the same time-dependence, while Kim{t,r) must 
have sin(fcf) time dependence. But, by Eq.(5.21), /i2£m must then have sm{kt) 
time-dependence. Similarly, from Eq.(5.22), given that h[^^^^ has sin(A:t) time- 
dependence, /ii7m must also have sin(fct) time-dependence, as must Aim{t,r) . 
These conclusions are consistent with our choice of boundary conditions (7.3,4). 
Noting Eqs. (5.37,42,43), the Dirichlet conditions (7.3,4) are equivalent to the 
boundary condition (7.1), which is analogous to a specification of momenta 
in a {^2im ' '^t^iim ) rt;presentation. This also accounts for the minus sign in 
Eq.(5.40). 

For large r , the potential term in Eq.(5.44) vanishes sufficiently rapidly that 

£,2em becomes a superposition of outgoing and ingoing waves at radial infinity. 

Note that Q^"'' also obeys Eq.(5.44); thus, Eq.(5.42) in the RW gauge tells us 

that h['^^^ = r Q^m^^^" — ^('') large r . Now, the field equation (5.27) 
implies that 

That is, odd-parity metric perturbations diverge at large r , in the RW gauge. 
This is only a coordinate effect, as the Riemann-curvature invariants decay 
at a rate 0{r~^) at large r [42,43]. (A similar phenomenon occurs for the 
even-parity perturbations in the RW gauge.) Here, in the odd-parity case, we 
construct a gauge transformation to an asymptotically-flat (AF) gauge, in which 
the radiative behaviour of the metric perturbations becomes manifest. 
Our odd-parity AF gauge is chosen such that 

hteL^^it,r) = . (7.9) 
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Thus, in terms of the preceding RW gauge: 



u{-)AF _ Q _ 
"'Oem — U — ll^g^^^ 



A-)AF 
''Urn 
A-)AF 
''2im 



h 



(-)mv 



lim 
2 A/m 



drA 



dtAim 

2A 



(7.10) 

(7.11) 
(7.12) 



Given h'^j,2f'^ and h^^^J^^ as a starting-point, one can, from the above, de- 



termine A^„(t,r) and hence h^^l^^ and hl^^J^^ ■ On substituting for h'l^J^ 
from Eq.(7.11) into Eq.(5.24), one finds 



(-)AF 



■(-)RW 



''Um 



2Ae ° (-)w 

„2 "■l^m 



(7.13) 



Now, following the approach used throughout when studying boundary condi- 
tions at the final surface Hp (t = T) , set: 



sm{kT) 



(7.14) 



where the {h[^l'^^{r)} are real radial functions. 

The property a^7m ~ i~^)"^(^^-ke -m ^olds for the coefficients. Similarly, 
one can construct a corresponding expansion for Q^^^^ {t ,r) . Then, from 
Eq.(7.13), one has 



2A 



(r) 



sin(fct) 



(7.15) 



p -^M V ; sin(fcr) ' 
which is 0(r~^) at large r , as required. On using Eq.(5.37), one further finds 



-f 



where 
and where 
satisfies 



^^l^a2Um^2M ^^^sin(fcT) 



A-)AF 
^2ke 



(r) = rh[-}^^{r) 



{-)AF' 

2ke 



(-) 



-)AF 



(7.16) 

(7.17) 
(7.18) 
(7.19) 



As in the spin-0 case [33,34] and in the spin-1 case above, we have, for fc > : 



i2ki-{r) - rjeikr) , 



as r 



(7.20) 
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^2ke-{r) ~ ((^2fe«-) exp(iA;r;) + (^Z2ke-) exp(-ifcr*)), as r* oo(7.21) 

where the je{z) arc spherical Bessel functions, and where r* is the Schwarzschild 
Regge- Wheeler coordinate [49,63] of Eq.(4.38). Thence, one deduces the nor- 
malisation property 

/■Roo 

/ dr e« Calf- (r) ^^k'e- (r) = 2n Wu- f (s{k , k')+5{k , -k')) . (7.22) 

Jo Sj? \ / 

The resulting form of the classical action for odd-parity (spin-2) gravitational 
perturbations can then be expressed as a functional of the complex quantities 
{aikim-} which encode the boundary data on J^p for the odd-parity gravita- 
tional perturbations. Here, 



^ilLs [{«2fe«m-}] 



1 (£ — 2)\ f°° 

= T^Yl Yl 77T^ / dk k \z2ke-\^ \{a2kem-) - {a2,-kem-)\^ cot{kT) . 
1^ e=2m=-e ^'^ + ^>-Jo 

(7.23) 

From this expression, one proceeds as in [33,34] (for spin-0) and as above 
(for spin-1) to study the semi- classical quantum amplitude or wave function, 
proportional to exp(iS'^^^^gg) , as a function of the complexified time- interval 
T = |T|exp(— for < < n/2. Just as in the spin-0 and spin-1 case, 
one straightforwardly recovers the complex Lorentzian amplitude for odd-parity 
gravitational perturbations, on taking the limit ^ — > 0+ . 



8 Regge- Wheeler formalism — even-parity grav- 
itational perturbations 

Working with even-parity gravitational perturbations in the RW formalism is 
notoriously more difficult than working with those of odd parity. Yet, Chan- 
drasekhar [57] showed that solutions to Zerilli's even-parity equation [64,70] 
[Eq.(9.10) below] can be expressed in terms of the odd-parity solutions. One 
might therefore expect that our results for the even-parity action should mirror 
those for the odd-parity action. 

We expand the even-parity perturbations as 



(h^A (x) - h[+lit , r) + Hiemit , r) e 



,(o-6)/2 



(/2)i 



r^Ke^{t,r) (/a) 



(8.1) 



Here, the non-zero components of the (un-normalised) basis of the symmetric 
tensor spherical harmonics [(/i,2,3,4)ii]^^ are defined by 



(8.2) 
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[{fl)r4>)hm = 


-■ {d4,Yem) 
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[if 3)00] 


= Yem , 






[ih)ee]g^ = 


{def Yim 




, - (cot d)d^ 



[iU)4>4>]em = [id4,f + [sine cos e)de]Yer, 
The non-zero inner products are given by 

26-" 
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I 



d^ [ihy']em [imim^ 

d^ Uh]e,nUy%m^ 
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-2a 



e 5u' Smm' ) 
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8.3) 
8.4) 
8.5) 
8.6) 

8.7) 
8.8) 
8.9) 



(8.10) 
(8.11) 
(8.12) 
(8.13) 
(8.14) 



where we define =£{£+!). The even-parity basis is also orthogonal to the 
odd-parity basis of Sec. 4. 

Further, for the even-parity perturbed shift, one can write 









tm 



(8.15) 



Bxtm^i ,T)Yim , hi-^l{t,r){deYem) , hi+l{t ,r){d^Yi„ 
For the perturbed lapse, 

[N^'^'^+^],^= -lHoUt,r)e-''/'Y,m ■ (8.16) 

Again, -ffofon = (^1)™ -f^o£.-m , etc. Hence, for the linear-order perturbation 
/i^^J of the 4-dimensional metric, the quantities h'"^ , h^) and h^^^ behave as 
scalars under rotations (their odd-parity part vanishes), while h^}g , , h^]a 



and /i^y behave as vectors. For 0,6 = 6,(j) , the 2x2 angular block h^^^ is a 
tensor under rotations. The even-parity gravitational momentum components 
can, correspondingly, be written in the form 



em 



+ 



PKim{t , r) [{f3)ij] + PaUt . r) [{h)ij\ 



5.17) 
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Again, one can easily show that the P's in Eq.(8.17) are related to hi , H2 , K 
and G of the corresponding Eq.(8.1) by: 




+ ^e{i+l)e-/^(dtG(m) . (8.21) 



For even-parity gravitational perturbations, gauge transformations are in- 
duced by even-parity gauge vector fields (^^"'"^'*)^m , of the form: 

(8.22) 

Within the adiabatic approximation, these induce the following even-parity 



gauge transformations: 

H;,,^ = Ho,m-a'x[tl + 2[dtX^+l) > (8-23) 

H[,^ = Hum - e-»(a.x(+l) _ e^[da[ti) > (8-24) 

H!,,^ = H2em- a'x[+l-2[drXi+l) , (8.25) 

K',^ = Kem - (^^) , (8.26) 

G'em = - 2-'^2fti ' (8.27) 

41 = 41+e-"^^li-r^(atXW) , (8.28) 

h[tl = h[Z-e^x[tl-r'{drXitl) ■ (8.29) 



As in the odd-parity case, we would like to construct gauge-invariant com- 



37 



binations of components of the perturbed 3-geometry. Following [68], we define 

kura = Kira + e-''(^r[drGern)-(^^^jj , (8.30) 

k2em = \[e''H2tm-e'''^dr(re'''^Kim)^ . (8.31) 

It can be shown that both the functions kx^m and k2em are gauge-invariant. 
For future use, in the calculation of the even-parity classical action, we define 
[73] the linear combination of kum and k2im '■ 

re{l + l)kum + 4re-2« k2em ■ (8.32) 

At this stage, as with the odd-parity case, we again make use of the property 
of the uniqueness of the (even-parity) RW gauge. In the RW gauge, one has 

^ V / (8.33) 

-2e''(5,4ti) , 

Then, 



ttRW _ tt I „2 



+ r(l + ^ra'^ (dtGem) - a' h] 



(+) 
oem 



5.34) 



Next, 



H^Z = Hum + [a' - e--^ (/I'+i - (a.G,„) 
+ r^e-4{^rfG,m-2^r{^^^\ . 



(8.35) 



Further, 



with 

together with 
and 



= K.ra - (^) (/llJi - y (drG.ra)) (8.36) 

GZ = = G,„-2XW , (8.37) 
4^:^^ = , (8.38) 



/^It:^^ = , (8.39) 
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where (in the RW gauge) 



4tl = ^-\^r\dtG,^)-hl,+ij , (8.40) 

^iti = e-''(h[+l-lr^[drGem)^ , (8.41) 

X^tl = \Gim . (8.42) 

At late times, following gravitational collapse to a black hole, in the ab- 
sence of background matter and in the adiabatic approximation, the even- 
parity RW field equations are seven coupled equations for the four unknowns 
{H^^,H^^,H^^,Kj^^). Assuming that £ > 2 — that is, that we are 
studying dynamical modes — we give here those RW field equations which are 
of first order in r and t [70] . These are, respectively, the {tO) , {tr) and (rB) 
components of the linearised field equations; 



(d.H^^)+f^e<^H^^ = e'^d,(K^;^ + H^Z) , (8-43) 



-re«(^l-^j(a,i^,T) , (8.44) 
= (a.ifoZ) - (drKZ') + ^-^H^Z , (8-45) 

and the (9^) component 

ttRW _ ttRW _ ttRW fd AR\ 

We also give one second-order equation, namely, the (rr) component: 



e / TO 

r \ r 



(8.47) 



Following Eq.(8.44), we find, for the gauge-invariant component defined in 
Eq.(8.32): 

{dtqun^) = t{t+l)(r(dtKZ)-e-''H^Z) ■ (8-48) 

We also find 

(Swwm) = - ^(^ + I)(2e-'^fc2te+ (l + ^ra'JfcH^I , (8.49) 
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where, in the RW gauge, 

qum = 2re-'^HP^ -2r^e-'^{drK^^)+2{Xr + 3m)KP^ , (8.50) 

with A = |(£ + 2)(£ — 1) . We can now solve for kum , k2tm in terms of qum 
and its radial derivative {drqum) , giving 

= 2(Ar + 3m) (( '^i+irO +^"-) ' ^^'''^ 
e°- I ( 3m \ \ 

^2^'" = - 4(Ar + 3m) H^^'^"'") + " ~ j ^""^ ) ' ^^'^^^ 

Further, Hf^ and Kf^ can also be written in terms of kum and kiim ■ 

9 Classical action and boundary conditions — 
even-parity gravitational perturbations 

As in the case of odd-parity gravitational perturbations, we can exploit the 

uniqueness of the RW gauge to simplify the even-parity action and to obtain a 
general gauge-invariant form for the even-parity classical action S\^!^^ . In the 
RW gauge, this is 



o(2) 
class 



"'^^ U'mm' 
1 \ " /"^°° ttRW* tjRW >o1^RV 



RW* uRW 
Ktm 



(9.1) 

Again, we would like to put the action into the form J drip (dtip) , where tp 
is gauge-invariant and obeys a decoupled wave equation. Since qum is the 
only unconstrained gauge-invariant even-parity quantity which involves only 
perturbations of the intrinsic 3-geometry, one might expect that Eq.(9.1) should 
reduce to the form 



o(2) 
class 



{qum} = ^E E / dr(numqum + {drZem))\ , (9.2) 



for some variable Zim , where TTUm is the gauge-invariant momentum conjugate 
to qiern ■ This is in fact the case. First, make in Eq.(9.1) the substitutions (as 
mentioned at the end of Sec. 8) for H^^^ and in terms of ku„i and k2im ] 

then substitute the expressions (8.51,52) for kum and k2im in terms of qum ■ 
After several integrations by parts, we arrive at an action of the form (9.2), with 

^ I pi?H' 1 3m\^jy ^3a/2l 



2(Ar + 3to) 



r 
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(Ar + 3rn) \^ (A + l) 

2(Ar + 3m) I (A + l) 



i^^Xe^/M , (9.3) 



7, - ^3 a/2 pW Tv-W , gl^" -^fen ^ pW a/2 Vq 



pW _ a/2(^pRW _apW pW ,c.r.\ 

This expression for the even-parity classical action simplifies yet further, 
since the linearised field equations imply that 



^-.r = ^e«ifC • (9.6) 



r 



Further, Eqs. (8.43,44) show, with the help of Eq.(8.48), that 

Are" 

= 2(Ar + 3m)^^*'^"-^ ' ^'"'^ 

('2') 

Eq.(9.2) for the even-parity S).-^!^^ then reduces to an expression of the desired 
form: 



C(2) 

class 



(9.8) 



where ^2^?! is defined as 

^2^™ - (Ar + 3m) ' 

We have made use of the assumption above that the specified perturbations 

ft-l^'' of the spatial 3-metric on the final boundary have been taken to be 
real. Of course, for the Dirichlet boundary-value problem with T rotated into 
the complex, the classical solution for the metric and scalar field will have both 
an imaginary part and a real part. 

Given the uniqueness of the RW gauge for even-parity modes, one can see 

(2) 

that Eq.(9.8) for S^^^^^ is in fact valid in any gauge, with a vanishing contribution 
from the total divergence since qum , as given by Eq.(8.28), and therefore also 
^2£m ' gauge- invariant. There are obvious similarities between Eq.(9.8) and 
the classical massless-scalar action of [33,34], with ^^'^ Co£m+ differing 

only by an ^-dependent normalisation factor. This should not be surprising, as 
scalar spherical harmonics have even parity. 

Again, one can show that the gauge-invariant quantity ^2eL satifies Zerilli's 
equation [70] 

e-''dr{e-'^{d4Z))-{dt)'ciZ-Vtkitl = , (9.10) 
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(+) ^ ( 1 _ ) ^ i. > . 

\ r J r3(Ar + 3m) 

(9.11) 

Now, both for odd and even parity, the field equations for the metric perturba- 
tions have been reduced to the two wave equations (5.44) and (9.10). 

In contrast to the odd-parity case, where we assumed an initially stationary 
multipole, here for even parity we treat ^2fei by analogy with the massless- 
scalar-field quantity ^oem+ , and impose the Dirichlet boundary condition 

dti(0,r-) = (9.12) 

at the initial surface S j (i = 0) . Proceeding now by analogy with the separation- 

of-variables analysis of Sec. 7 for the odd-parity case, we find that, if has 
s,m.{kt) time-dependence, then so must also, whereas -f^^m i^ust have 

cos(fcf) time-dependence. Consistency with the gauge transformations (8.33- 
39) implies that these time dependences are valid in an arbitrary gauge, and 
further that G(,rn and h^-^^ have sin(fct) time dependence, whereas /iq^^ has 
cos(fcf) time dependence. Consequently, qum must have sin(fct) time de- 
pendence, whence the boundary condition (9.12) is justified through Eq.(9.9). 
(Alternatively, one could instead have studied normal- mode time dependence.) 
Following the scalar-field analysis of [33,34], we can write 

_ r ^uJ+) sin(fct) 

sin(fcT) 



2m 



\ 2(^A2(A + I);-'' + .'U-//);-2 + 9A//)2,. + ()„,3^ 



e-iit,r) = *4:L^^:i(r)^:^ , (9.13) 




where the {02^]^} are suitable even-parity 'Fourier coefficients', and where 
{^2fejm(^)} ^^^^ radial functions. These functions satisfy 

e-''^^ + {k'-Vt\r)y^J = . (9.14) 

Regularity at the origin implies that 

^iVeir) - (const.) X rje{kr) (9.15) 
for small r . Again, at large r , the potential vanishes sufficiently rapidly that 

^2ke asymptotic form 

^itkr) ~ ((4+i)exp(ifcr:)+(4+)*)exp(-ifcr:)) , (9.16) 

where {^2^}^} complex constants. Then, the classical action S^^^^ for even- 
parity gravitational perturbations reads 

c(2) rr„(+) il 

■^class [t"2/s£mJ J 

(9.17) 
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where the notation is in hne with that for spin-0 and for odd-parity fields. The 
coordinates {a2kem+} label the configuration in fc-space of the even-parity part 
of the metric perturbations on the final surface S p . 

Let us now rc-asscmblc both the odd-and even-parity metric perturbations. 
As above, we consider for simplicity odd-parity metric perturbations which are 
initially static (Neumann problem) and even-parity metric perturbations which 
vanish initially (Dirichlet problem), on the space- like hypersurface S/ . The 
total classical spin-2 action is then 



oo e 



^TfiX] X] H (f + 2V / ^^^l^2fc£^l^l("2fc£™^) + (-Pa2,-fe«mp)|^cot(A;T) , 

(9.18) 

where the complex coefficients {a2Mmp} obey 

a2kemP = -P (-1)™ tt2,-fcf,-mP • (9.19) 

Here, P takes the value ±1 according as the parity is even or odd. 

As in the case of odd-parity metric perturbations (Sec. 7), the even-parity 
metric pertiirbations also diverge at large r , except in a special gauge, the 
asymptotically-flat (AF) gauge. In the AF gauge for even parity, as for odd 

parity, all physical components ^[^^(i,) = \'y^'^'y'"^\^ h\ju (that is, all compo- 
nents of projected onto the legs of a pseudo-orthonormal tetrad oriented 
along the unperturbed it ,r ,9 ,(j)) directions) fall off in the wave zone more 
rapidly than r^^ , except for the transverse (angular) components, which carry 
information about the gravitational radiation. In the new (AF) gauge, for even 
parity, one has 

= = HtL = . (9.20) 

Further, from the even-parity gauge transformations (8.23-29), one has 

= H^Z-a'x[^l + 2(d,X^+l) , (9.21) 



= H^Z + e-''{drXi,ti)-^^{dtXiti) , (9.22) 
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ttRW 
^Um 
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-'-(drX, 
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Kt^ = i^ir - ( =^ ) , (9.24) 



{drXitl) , (9.23) 



'-'tm — ~ 2A^2fct ' (9.25) 



= e-'^4ti-r'{dtX^tl) , (9.26) 
h'.'elf' = -e'^ntl-r'{drXitl) , (9.27) 
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where a hat denotes a gauge function in the AF gauge. Therefore, once given 
^tof . HiZ and K^^ , then Eqs.(9.21,22) can be solved for Jt^+l and x[+l . 
Thence, Eq.(9.26) can be used, in order to solve for ^2ftn • solving these 
equations, one chooses the arbitrary functions which arise such that asymptotic 
flatness is still satisfled. Thus, the AF gauge is consistent. 



10 Conclusion 

In this paper, we have taken over the scalar (spin-0) calculations of [33,34], with 
the help of the angular harmonics of Regge and Wheeler [63], to include the 
more complicated Maxwell (spin-1) and linearised graviton (spin-2) cases. For 
spin-1, the linearised Maxwell field splits into a part with even parity and a part 
with odd parity; a different treatment is needed for each of these two cases. In 
both cases the relevant boundary conditions involve fixing the magnetic field on 
the initial space-like boundary S/ and final boundary Tip . The main result is 
an explicit expression (6.19) for the classical (linearised) Maxwell action, as a 
functional of the final magnetic field, subject to the simplifying assumption that 
the magnetic field on the initial surface S/ is zero. From this, the Lorentzian 
quantum amplitude for photon final data can be derived, as in [34] for spin-0 
perturbative final data, by taking the limit 9 ^ 0^ oi exp(i5'ciass) ) where S'dass 
is the action of the classical solution of the boundary-value problem with pre- 
scribed initial and final data, with complexified time- interval T = \T\ exp(— i9) , 
where <9 <tt/2. 

Linearised gravitational-wave (s = 2) perturbations about a spherically- 
symmetric Einstein/massless-scalar collapse to a black hole have also been stud- 
ied here. As for Maxwell (s = 1) perturbations, the principal aims for s = 2 also 
are (1) to specify suitable perturbative boundary data on the final space-like 
hypersurface Tp at a late time T, subject (for simplicity) to the initial bound- 
ary data on S/ (time t — 0) being spherically symmetric; (2) to express the 
spin-2 Lorentzian classical action S'dass as an explicit functional of the 'suitable' 
boundary data above, and of the proper-time interval T , once T has been ro- 
tated into the complex: T — > |T| exp(— i9) , for < 9 < 7r/2 ; (3) given S'dass , 
to compute, following Feynman, the quantum amplitude for the weak-field final 
data, by taking the limit of the semi-classical amplitude (const.) x exp(iSdass) 
as ^ 0+ . 

As in the s = 1 case, it is also necessary for s = 2 to decompose the metric 
perturbations into parts with odd and even parity. The main difference on 
moving from the s = 1 to the s = 2 case is a considerable increase in algebraic 
or analytic complexity, to be expected since one deals with tensor fields rather 
than vector fields. 

Some indications towards unification of these ideas for perturbative fields 
of different spin s appear already in Sees. 2,3. For s = 1 , the quantity most 
naturally specified as an argument of the quantum wave-functional, on a bound- 
ing hypersurface such as E^? , is the (spatial) magnetic field Bi , subject to the 
condition ^Vfe.B*' = . Correspondingly, for linearised gravitational waves 
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(s — 2), the natural boundary data were found to be the (symmetric, trace- 
free) magnetic part Hi^ of the Weyl tensor [52,53], subject to ^SI = . In 
2-component spinor language, these correspond (s = 1) to a particular 'projec- 
tion' of the (complex) symmetric Maxwell field-strength spinor (^ab = 4'{ab) j 
and (s = 2) to a corresponding projection of the totally-symmetric (complex) 
Weyl spinor "^abcd = ^(abcd) ■ Of course, as treated in [35], these boundary 
conditions constructed from (pAB and abcd are special cases of the natural 
boundary conditions for gauged supergravity [71-73]. In our work on quan- 
tum amplitudes for spin-i [35], the natural boundary conditions also involved 
a corresponding projection of the spin-i field. Although 2-component spinor 
language might (to some people) seem a luxury in treating bosonic fields de- 
scribing photons or gravitons, above, it is practically a necessity in treating the 
corresponding fermionic (massless) neutrino spin-i field, as in [35], and (for 
supergravity) the gravitino spin-| field, on which work is in progress [36]. 
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